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DISCUSSION OF ANALYSIS OF GROUND-WATER LOWERING 
ADJACENT TO OPEN WATER 
PROCEEDINGS-SEPARATE 106 


Marruew I. RorapauGu,t A. M. ASCE.—An orderly analysis of ground- 
water lowering adjacent to open water is presented in this paper, the entire 
analysis being based on the assumption that the total discharge is uniformly 
distributed among an infinite number of wells in a ring or, in other words, 
that each wellpoint is pumped at the same rate. The theoretical water level 
within the ring was found to dip away from the river; experimental data were 
presented indicating a dip toward the river. This difference was explained 
as probably being due to the pump location and pipe friction. To the writer 
it seems that the analysis is based on an assumption that does not fit the 
problem and that the difference between theoretical and observed results is 
only partly the result of pipe friction. 
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In practice, a suction pump is connected to a line or ring of wellpoints. 
The system is one of constant drawdown at all points on the ring, not uni- 
formly distributed pumping. An analysis based on constant drawdown will 
be consistent with the field installation and will produce a solution that fits 
the problem—that is, it is presumed that the excavation is level and that 
a horizontal, rather than a dipping, water level is desired. 

Fig. 15 shows a flow net for the case of a pumping well of very small dia- 
meter located near a surface source of water. The flow lines are ares of circles 
having their centers on the line source. The lines of equal drawdown are 
nonconcentrie circles having their centers on a line normal to the line source. 


4 Dist. Engr., Ground Water Branch, U. 8. Geological Survey, Louisville, Ky. 
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The theory and equations for this system (for the artesian case) have been 
published by Morris Muskat® and C. E. Jacob.® 

This system is readily adaptable to the problem of a ring pumped at con 
stant drawdown. All that must be done is to place the ring on one of the 
circles of equal drawdown rather than concentric with the hydraulie center 
of the system. Pumping from the ring will not be uniformly distributed, 
being greater on the riverward side where gradients are steeper and stream 
lines are more dense, and less on the landward side where gradients are less 
steep and stream lines are less dense. The flow pattern outside the ring will 
be the same for pumping the ring as for pumping a small well at the same 
rate, if the small well is located at the hydraulic center. 

The drawdown at any point in the field for the case of a small-radius 
pumped well is expressed by the equation, 


in which qg is the pumping rate; a, the distance from the line source to the 
center of the well; and r, the distance from the pumped well (or hydraulic 
center of flow system) to any point at which drawdown is desired. Other 
terms are as defined in the paper and as shown in Fig. 16. For the riverward 
profile, 8 = 0, and 


q 


For the landward profile, @ = 180°, and 


q 
awk B 


Eqs. 23 and 24 are similar to Eqs. lla and lle except that the former are 
for a single small-diameter pumped well located a distance a from the source, 
whereas the latter are for a ring being pumped at a uniformly distributed dis- 
charge, the center of the ring being at a distance p from the source. 

The location of the center of the ring with respect to the hydraulic center 
of the flow net is determined as follows: Let the center of the ring be at a 
distance p from the source. Then, the eccentricity. e, is p — a. 

The term e can be evaluated by writing Eqs. 23 and 24 for two points on 
the specified profiles at locations diametrically opposite each other on a ring 
of radius r;; thus, in Eq. 23, d = r, — e; in Eq. 24, d’ =r, +e. Since the 
drawdown is the same at all points on the ring, the log terms of these two 
equations are equal. From this relationship, 


Since e = p — a, then , 


*“The Flow of Homogeneous Fluids Through Porous Media,"’ by Morris Muskat, McGraw-Hill 
Book Co., Inc., New York, N. Y., 1937 p. 175. 

“Engineering Hydraulies,"’ Ed. by Hunter Rouse (‘'Flow of Ground Water,"’ by C. E. Jacob, Chap- 
ter 5), John Wiley & Sons., Inc., New York, N. Y., 1950, p. 345. 
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Thus, if the location of the center of the ring, p, and the radius of the ring, ri, 
ure known, the location of the hydraulic center a is determined. Knowing 
a, the profiles on the riverward and landward sides are computed from Eqs. 


23 and 24. 
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The equation for discharge for a given drawdown at any point on the 
ring or any point within the ring is determined from Eq. 23 or Eq. 24 by 
substituting d =r, —e; =r, e=p—Vp—ry; a= 
Za = 24 = hg = he’ = andz, = H — h,; thus— 


wrk(2Hh, — h?, ‘ 
pt+tvp-r, 
log 


Eq. 27 is similar to Eq. 6 except that eccentricity is included. For a single 
well or a very small ring, when r; is small relative to p, the log term approaches 


2 P as in Eq. 6. 
ry 
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The significance of the constant drawdown approach is demonstrated in 
Fig. 17, which shows profile A, Fig. 11, and the solution for the same case by 
the constant drawdown method. Note that the constant drawdown solution 
produces a horizontal profile within the ring which is consistent with the results 
desired and is consistent with the field condition of a group of wellpoints 
connected to the same pump. It should be emphasized that, for this partic- 
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Fic. 17.--Drawvown Curves sy Two Metnops 


ular case, the eccentricity is small, so that the discharge required to produce 
the same center drawdown, as determined by Eq. 27, is only about 1.0% 
larger than that determined from Eq. 6. 

The discharge as computed by the two methods differs only in the log term. 
The discharge required to produce the same center drawdown for the constant 
drawdown method will be larger than that by the uniformly distributed 


discharge method by the ratio —----—--————._ As r; approaches p (that 
ry 

is, for excavations very close to the surface source) this ratio will become large; 
for p = r, the ring becomes tangent to the source, and both methods break 
down since this condition would infer direct pumping from the surface source. 

The differences in the profiles outside the ring will become larger as r; 
approaches p, as shown by comparing Eqs. Ila and Ile with Eqs. 23 and 24. 
For the constant drawdown method, as the ratio r,;/p is made larger the ec- 
centricity is increased, a becomes smaller (the hydraulic center moves river- 
ward), and the riverward profile becomes steeper. The profiles outside the ring 
based on Eqs. lla and Ile are independent of r; so that changing r; does not 
alter the profiles. 
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This discussion has been limited to the problem of a single ring; but the 
criticism of the basic assumption of a ring having equal distribution of dis- 
charge rather than a ring of constant drawdown applies to all the variations 
covered in the paper. 

Two other points might be mentioned briefly: (1) The location of the line 
source and (2) the effects of river-water temperature. 

1. Location of the Line Source.—In problems of water-supply develop- 
ment based on induced infiltration of river water at locations along the Ohio 
River, it has been determined that the line source is not at the shore line but 
is usually several hundred feet offshore. The theory assumes a vertical 
boundary between the source and the aquifer. In nature, water is entering 
the aquifer vertically in the river bed and then flowing horizontally toward 
the well. The problem is handled by replacing the complicated flow condi- 
tions under the river by an “effective’’ line source located at a position off 
shore so that the flow pattern in the well field is approximated. In the case 
of water-supply development, location of the ‘effective’ line source is deter- 
mined by controlled field pumping tests, inasmuch as its location may be im- 
portant in establishing the minimum yield of the installation. Assumption of 
the shore line as the line source in problems of unwatering is on the safe side; 
that is, the actual pumping rate should prove to be less than computed. 

Effects of River-Water Temperature.—The author uses the definition: “k is 
the effective coefficient of permeability.” Normally ground-water tempera- 
ture is nearly constant and the viscosity does not vary greatly, so that small 
errors are lumped into the permeability term. The viscosity of water in the 
range encountered in the problem in question varies about 1.5% per degree 
Fahrenheit; thus a rise of 1° F will increase flow rates about 1.5%. The 
temperature of the surface source becomes an important item in the problem 
if unwatering is continued over a long time. As an example of the range 
involved, Ohio River water varies in temperature from 32° to about 85° during 
the year. Continued pumping near the river induces water of varying tem- 
perature to enter the aquifer and flow to the point of pumping. The problem 
becomes complicated as river water mixes with ground water of a different 
temperature, heat exchange occurs between the water and sand, and flow 
rates at each point adjust themselves to the temperature prevailing at that 
point. The distance from the river is an important item. Detailed records 
at an installation at Louisville, Ky., for the period from 1945 to 1951 show 
a range in temperature from 47° to 64° for the discharged water. Statistical 
analysis shows the flow rate for a condition of constant river level and con- 
stant pumping level to vary from 85% to 115% of the average value. In 
the case of unwatering over a long period, the pumping rate might vary as 
much as 30% or 40% between summer and winter. Also, this factor will 
cause variations in the pumping distribution in the ring—in the summer, 
discharge from wellpoints on the riverward side will increase and in the winter, 
discharge on the riverward side will decrease. Perhaps Mr. Avery can furnish 
information as to how long pumping continued in his example, what seasons 
it covered, whether temperatures of discharged water were recorded, and 
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whether any differences in pumping rate or water level were noted which 
would be related to the temperature variable. 

Inasmuch as the temperature variable was neglected in the analysis, the 
six assumptions given at the beginning of the paper should be increased to 
seven by adding (g) The river temperature and groundwater temperature 
are equal and do not change during the pumping period. 

The conclusion that errors as great as 10% to 15% are to be expected 
might be questioned. Where pumping tests have been run to evaluate the 
“effective’’ distance to the line source and the permeability, and where tem- 
perature has been included in the computations, water supplies based on in- 
duced infiltration generally have been predicted with this degree of accuracy. 
However, if the line source is assumed as the shore line and the temperature is 
neglected, greater errors should be expected. 

Summary.—On the basis of the assumptions made, the paper presents 
a sound, orderly treatment of the problem. 

The assumption that discharge is uniformly distributed in the ring does 
not fit the problem, inasmuch as a sloping profile is produced by the theory 
whereas a horizontal profile is desired. Analysis based on the assumption 
of constant drawdown at the ring will produce the desired horizontal profile 
and will be consistent with the field operation of pumping at constant draw- 
down. 

Observed results confirm the view that the constant drawdown approach 
gives a closer check between theory and practice. 

I:ffects of river temperature were not considered in the analysis. Neg- 
lecting this item may introduce errors of as much as 30% or 40%. This 
item alone is considerably larger than the 10% to 15% limit of error assumed 
by the author. 


8. J. Jounson,’? A. M. ASCE.—The dearth of information on the engineer- 
ing aspects of ground-water lowering for construction purposes and the need 
for such information combine to make this paper one of unusual timeliness and 
importance. In the United States, ground-water lowering for construction 
purposes has been all too often resorted to as an expediency when in trouble 
or used without recognition of its possibilities and limitations. The practical 
aspects of ground-water lowering operations are of primary importance and an 
impressive quantity of working knowledge is available. However, the the- 
oretical aspects have not been developed simultaneously with the practical 
aspects, even though theoretical considerations can supplement field experience 
to advantage and can also form a logical framework for the acquisition and 
organization of practical knowledge and experience. Publication of this paper, 
reviewing some of the available theoretical methods, is therefore of significant 
value. It is unfortunate that the scope could not have been extended to in- 
clude a complete review of ground-water lowering theory. 

It would have been of some interest and importance if the author had re- 
viewed the status of the available knowledge on ground-water lowering methods 


'Chf., Embankment and Foundation Branch, Waterways Experiment Station, Corps of Engrs., 
U.S. Dept. of the Army, Vicksburg, Miss. 
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of analysis. The writer’s understanding of this subject is as follows: The 
equation for gravity flow to a single well, which depends on the validity of the 
Darcey law, was formulated by Mr. Dupuit in 1863,° as was stated by the 
author. Before 1900. Mr. Forchheimer developed the equations for flow to 
groups of wells and used the method of images to develop equations for the 
flow to wells supplied by a line source.*:!° The formulas for practically all cases 
of the flow of water to either well or wellpoint systems (either single or multiple 
stage with the wells arranged in a line, rectangle, or cirele) were presented prior 
to 1930. For example, the first edition of the book by W. Kyrieleis" on ground- 
water lowering by pumping from wells was published in 1913 and presented 
most of the equations needed for caleulating the effect of pumping from wells. 
The second edition, by Mr. Kyrieleis and W. Siehardt, was published in 1930. 
Other references are those by H. Weber and Mr. Sichardt"; the latter contains 
formulas for multiple-stage well systems supplied by a circular souree. Un- 
fortunately, practically none of the available theory is found in’ American 
engineering literature. Additional investigations are needed to develop further 
certain theoretical considerations, but the available theory is sufficient for many 
practical applications. Therefore, the significance of this paper is readily 
apparent. 

The Dupuit assumptions have been subjected to much criticism and many 
have questioned the usefulness of the theoretical methods for estimating draw- 
down which are based upon it, although its usefulness for satisfactorily estimat- 
ing the quantity of flow has been acknowledged. Mr. Casagrande has pointed 
out that the Dupuit assumptions are satisfaetory if the slope of the lowered 
water surface is not too great. This has been confirmed through model tests 
by H. E. Babbitt, M. ASCE, and D. H. Caldwell,'4 A. M. ASCE, in relaxation 
solutions by Shih-Te Yang," and by the model tests and relaxation studies of 
H. P. Hall.® Although extensive data are not available, it appears that the 
Dupuit assumptions do not involve an excessive error in the estimated draw- 
down at a distance beyond from ten to fifteen well diameters from the center of 
the well. Consequently, it can be assumed that the Dupuit-Forehheimer 
theory for steady-state gravity flow to well systems is valid for ordinary re- 
quirements in estimating the quantity of flow and the lowered ground-water 
level, except that the computed values of the loweved ground-water level near 


Les fontaines publiques de la Ville de Dijon,” by H. Darey, Dijon, France, 1856. 
die von Brunnenanlagen und Sichersehlitzen,”’ by Ph. Forchheimer, Zettechrift 
des und Vereins 2u Hannover, Vol. 32, 1885 
be: by Ph. Forchheimer, Zettschrift des oesterrlichiachen 
Inyenteur und Vol. 50, 1898. 
erabsenkung by W. Kyrieleis, Julius Springer, Berlin, 
Germany, 1928 
Phe Rochwete von Grand wasserabsenkungen mittles Rohrbrunnen,”’ by H. Weber, Julius Springer, 
Serlin, Germany, 1028 
Pas Passungsvermogen von Rohrbrunnen und seine Bedeutung fir die Grundwasserabsenkung 
insbesondere fur grossere \bsenkungstiefen,”” by W. Siehardt, Julius Springer, Berlin, Germany, 1928 
Pree Surface Around, and Interference Between, (iravity Wells,” by H. bk. Babbitt and D. H. 
Caldwell, Bulletin No. 30, Univ. of Hiinois, Urbana, I1., Jan 7, 1948 
scopage Toward a Well Analyzed by the Relaxation Method," by Shih-Te Yang, thesis presented 
to Harvard University, in Cambridge, Mass., in 1949, in partial fulfilment of the requirements for the degree 
of Doetor of Philosophy 
1 Tavestigation of Steady Flow Towards a Gravity Well,” by H. P. Hall, thesis presented to Harvard 
University, in Cambridge, Mass., in 1950, in partial fulfilment of the requirements for the degree of Doctor 
of Philosophy 
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the wells will be considerably too low. However, this does not affect the prac- 
tical usefulness of the theory, since it is the lowering at a distance from the 
wells—as at the center of an excavation—that is the practical interest. 

The author's procedure in developing Eq. 3 is not clear and the resulting 
effect is that his mathematics appear to be incorrect. Apparently, this is be- 
cause of his desire to conserve space, but the reader cannot readily supply the 
missing steps. The author assumes for simplicity that there are an infinite 
number of wells in Eq. 3 and loses the reader in proceeding to Eq. 4. For 
equal well discharges, Eq. | becomes: 


if = b> log. S; — > loge Ra] (28) 


in which ¢ = 1,2,---n. For wells of diameter r; d@ so close that they touch, 


the number of wells, n, equals 2 4 r_orn= Ad so that the flow per well be- 
do 


comes oe , in which Q is the total flow from the system. Thus, introducing 


the expression for S, = y 4p? + — 4 pri cos 6, Eq. 28 becomes 


Q dé d0 
H? — 2, = x [| #r loge Si — log. 


An evaluation of > log, S, d@ is the integration referred to in Eq. 4, which re- 
sults in Iq. 5. 

The writer has found it convenient to express the terms in the brackets in 
eq. 28 in the forms of dimensionless expressions and graphs for a finite number 
of wells. If Q is the total well flow, Eq. 28 may be expressed as 


S 
H? — 2’, = Q [= log, 


t=n 
which becomes, for a circular well arrangement, since >> log, R, = n log, ry 


I? - 2, = log. Si loge |. 


i=l 


Using the notations from Fig. 4, 


n 


or . . (32) 


2 
4 (2) cos 1) | 
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log. Si = loge ri + loge ( P y -4 cos (i — 1) as | . . (33) 
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i=n 
The summation = > log. S; now can be expressed in dimensionless form 
i=l 


and plots of it can therefore be simply prepared. Eq. 31 thus becomes 
» 9 Q 
2, = — [ log, ri +L — log, ri] 
wk 
QOL 
r-2,= 
ak 


where 


2 
L= ¥ ton. [1 +4 (2) 4 F cos 1) 


The value of L is plotted in Fig. 18 for various values of p/r; and different num- 
bers of wells. The flow per well is assumed to be the same for all wells. It 
ean be seen that L is only slightly dependent on the number of wells used. 
Therefore, the expression, L = log, 2 p/ri, as given (in effect) by the author in 
Eq. 5, may be used regardless of the number of wells when p/r; is greater than 
about 2, and for a moderate number of wells when p/r; is less than 2. 

The equation for the lowered water surface in the wells, when a limited 
number of wells is used, can be obtained simply and is of interest. If the dis- 
tance from the impervious layer to the water surface in the well is denoted by 
2», it is apparent from Eq. | that for a circular arrangement of wells with an 
equal flow per well 


i=l 


i=<=n 

It is apparent that - > log, S, can be obtained from the previous evalu- 
i=l 

ation of this summation by simply letting 2 p equal the distance from the well 
at which the drawdown is desired to the center of the imaginary well system, 
which distance may be called 2p’. For most well systems the value of the 

summation can be approximated by log, 7: + log. —E Mr. Sichardt explains 

1 


that 


r 


1 1 1 
log. Ry = log, ry — — log. —— 


in which r, is the well radius."’ It follows that the evaluation of the water in 
the well may be found from the equation, 


9 
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The author limits his discussion to circular well arrangements. If a rec- 
tangular system of wells or wellpoints is used, it may be converted into an 
equivalent circular system having an enclosed area equal to the area of the 
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rectangle, or the rectangular system can be analyzed by using the graph pre- 
pared by Mr. Weber.”* The writer has found it convenient also to prepare dimen- 
sionless graphs of © log, S, and So log. R, for lines of wells. These graphs also 
can be used to find the lowering at any point when using a rectangular area 


Values of L, Eq. 35 
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Ratio, + 
Fig. 18. -Cinmet car wirn Link Source 


The example given by the author is particularly interesting beeause it 
demonstrates an engineering approach to a ground-water lowering operation. 
Unfortunately, there are only a limited number of cases where such an approach 
is actually used. The author did not indicate how he treated the fact that the 
wells only partly penetrate the previous stratum. This is a practical feature of 
considerable importance. He apparently assumed that the wells completely 
penetrated the pervious stratum, which is a different assumption from that 
sometimes made. The author’s viewpoint on this phase would be appreciated. 
The length and spacing of the wellpoints are pertinent and should, perhaps, be 
indicated. 
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In conclusion, the writer wishes to emphasize his belief in the substantial 
significance of this paper and in the practicability of applying this approach. 
At present, only a few engineers have been applying in practice the available 
theory for ground-water lowering. It is hoped that the paper will stimulate 
interest in this field. 


Howarp P. Hatt,'’ A. M. ASCE.—Convenient expressions are developed 
in this paper for determining the drawdown surface in a number of multiple- 
well problems that have not been treated in such detail before. The writer 
would like to outline an analytical approach that may simplify the derivation 
of the equations to some extent, and thereby facilitate further investigation. 

Fig. 19 illustrates the case of a single gravity well near open water in 
the form in which it was originally analyzed by Mr. Forchheimer.'® The 
ground-water surface before the beginning of operation of the well is assumed 
to have been horizontal at the elevation of the open water surface, the pervious 
stratum is assumed to be homogeneous and isotropic, and the well is assumed 
to penetrate the full depth of pervious material to an underlying horizontal 


SECTION 


Fig. 19. Forcuuemer SoLtution TO A SINGLE Near A River 


impervious base. Mr. Forchheimer imagined the body of open water to be 
replaced by a mirror image of the given well and the semi-infinite soil mass, 
and required that the imaginary well feed water into the pervious stratum at 
the saine rate at which it was being drawn out of it through the real well. 
Conditions of flow in the given soil mass would not be affected by this sub- 
stitution, but the problem could now be treated as a very simple well group. 
Making use of an analysis which he had already developed for a group of wells,'® 
Mr. Forchheimer derived the equation: 


id q 
2 = — 
us 


in which z is the elevation of an arbitrary point on the free surface, measured 
from the impervious base; // is the elevation of original ground-water surface, 
measured from the impervious base; g is the rate of flow from the imaginary 
well to the given well; & is the coefficient of permeability of the pervious 
stratum; and z, y are the distances from the arbitrary point to the imaginary 
well and given well, respectively. 

Referring now to the problem in the paper, let Fig. 20 represent a typical 
arrangement of a number of wells in a circle having a radius 7, and a center 
at a distance p from open water. Eq. 4 and Eq. 11 have been derived with 


17 Asst. Prof. of Civ. l-ng., Northwestern Univ., Evanston, Il 


18 “Grundwasserspiegel bei Brunnenaniagen,”’ by Ph. Forehheimer, Zettschrift des oesterreichischen 
Ingenieur und Architekten-Vereins, Vol. 50, 1898, pp. 629-648 


“Uber die Ergiebigkeit von Brunnenanlagen und Sicherschlitzen,"’ by Ph. Forchheimer, Zettschrift 
des Architekten und Inyenieuren Vereins zu Hannover, Vol. 32, 1886, pp. 539-564. 
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such an arrangement in mind, but with the well group idealized as a ring 
composed of an infinite number of small wells—that is, as a well having an 
annular cross section. Examination of the form of these equations suggests 
that the problem might be analyzed by replacing the given ring of wells with 


Fig. 20.— Ring of WELLS AND AN EQuivaALeNntT SINGLE WELL 


a single well having the same center and radius as the ring, and discharging 
at a rate equal to the sum of the rates of discharge of the given wells. Thus, 
if n is the number of given wells and q is the rate of discharge from each well, 
Eq. 39, applied to this substitute well, becomes— 


—for any point, M, outside the ring. 

Within the ring, Eq. 39 may also be used, but it must be noted that, 
theoretically, the drawdown contours within this area are concentric circular 
arcs whose venter is the center of the image ring. (The original analysis 
superposed the drawdown of the given well upon the negative drawdown of the 
image well.) Thus, the elevation of the drawdown surface at any point within 
the ring is the same as that of a point on the ring circumference at the same 
distance from the image center. Eq. 39 for a point, M’ (Fig. 20) within the 
ring may therefore be written: 


A comparison of Eqs. 40 and 41 with Eqs. 11 shows that they are the same. 
It should be noted that the introduction of the idea of an equivalent single 
well does not reduce the time and labor required for computing numerical 
results since the equations are the same, but the derivation is somewhat 
simpler since the analytical procedure is reduced to the direct application of 
one equation. 

The foregoing discussion indicates only that the equivalent single well 
may be substituted for a well of annular cross section without loss of theoretical 
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accuracy. It is also of importance to know to what extent.the number of wells 
ean be reduced and the radius of ring increased in relation to p (Fig. 20) 
before the device loses its theoretical validity. Fig. 11 contains drawdown 
curves for a ring of six wells, which can be compared to the curve for the 
idealized annular well, and, therefore, to the equivalent single well, under 
similar conditions. The differences indicated are not large. The following 
considerations extend this comparison a little further. 

Consider a group of n equal wells, each discharging at a rate qg, located 
on a circle whose radius is 7; and whose center is at a distance p from the shore 
line, as in Fig. 20. Let r, be the radius of a single well having the same center 
as the given ring and the same total rate of discharge as the n given wells, and 
producing the same drawdown at the center; that is, let r, be the radius of 
the equivalent single well corresponding to a ring composed of a finite number 
of wells. 

Beginning with the direct analysis, the following is obtained from Eq. 1: 


in which z, is the elevation of the free surface at the center of the ring, and 
z, is the distance from the center of the ring to the image of the nth well. 
An analysis of the same case by the equivalent single well gives 


e 


If the drawdown at the center is to be the same, it follows that 


le 
or : 
2P Te, °° Za 
(21 ( ) 
and 
= 
V Zi, » 


Eq. 46 is a convenient form in which to compare the radius of the equivalent 
single well to the radius of the given ring of wells. If p is large in relation 
to r, it follows that 2, = rz. = --- = 2, = 2p, and the ratio reduces to unity, 
as would be expected. As 7; approaches the magnitude of p, the corresponding 
changes in the ratio are insignificant until the circumference of the ring of 
wells comes very close to the shore line. For example, in the case of a ring 
of as few as four wells located at the quarter points, the radius of the equivalent 
single well is less than 1°, different from the ring radius until p becomes less 
than 1.1257,. Consequently, it appears that, from the point of view of 
drawdown at the center, the equivalent single well is an acceptable analytical 
device in most causes. 

The differences at the circumference of the ring are the greatest, as Fig. 11 
indicates. Lg. | applied to a point on the circumference of the ring may be 
written 


2. = q or Tn 7 
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The corresponding expression using the equivalent single well becomes 


n ¢ 
log. 


In cases where the drawdown is assumed to be not greater than 0.2 //, Mr. 
Casagrande recommends the following approximation as introducing an error 


of not more than 10%: 


« 
CASE II (EIGHT WELLS) At 


a’. 
Fic. 21.—Rario 2 Verses av Poin’ A 
mn ha 


in Which A is the drawdown = IT — 2. If Bq. 476 is divided by Eq. 47a with 
this approximation introduced, the resulting equation expresses the ratio of 


drawdowns computed by the two methods: 


Yr, °° Yo 


in which A, is the drawdown determined by use of the equivalent single well 
and hy is the drawdown determined by summation of the contributions of 
individual wells. 

iq. 48 provides the basis for the curves of Fig. 21. A ring of four wells 
and one of eight are considered separately. Each group is arranged so that 
the point on the ring circumference which is nearest the shore line is also mid- 
way between two wells. Thus, point A corresponds to the largest discrepancy 
between the equivalent single well solution and the conventional summation 
for a given number and location of wells. In view of the fact that these 
differences are localized, as the six-well drawdown curve of Fig. 11 indicates, 
it appears from Fig. 21 that, even at the ring circumference, the equivalent 
single well is rensonably reliable in most cases. 
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DISCUSSION OF REDUCTION IN SOIL STRENGTH 
WITH INCREASE IN DENSITY 
PROCEEDINGS-SEPARATE 228 


E. S. Barber,2 A. M. ASCE.-It is reasonable to expect pore pressures for 
compaction beyond the line of optimums since the phenomenon of an optimum 
moisture is caused by pore pressures developed in excess water which can- 
not move away quickly enough. Conversely, no optimum is obtained for very 
permeable soils or indeed for any soil under static loading if sufficient op- 
portunity for drainage is provided. 

Figure A corroborates the reduction in strength from overcompaction by 
impact; however, the continued increase in strength withdensity for static 
compaction indicates that something besides pore pressure is involved. Since 
shearing is restricted in static compaction, it is suggested that the shearing 
involved in impact may cause reduced strengths. 

To eliminate pore pressure as a factor, tests were made on dry sand. For 
a fine sand, with a 10-lb, hammer dropping 18 inches, a CBR of 5 as com- 
pacted was obtained while compaction to the same density by vibration gave 
a CBR of 18. Similar tests on a coarse sand gave CBR’s of 6 and 16, respec- 
tively. Apparently, the reduced density along the shear surfaces from impact 
is a major factor in reducing strength. 

A clayey sand and gravel, initially compacted by impact to 129 p.c.f. with 
8.8 percent water was subjected to 100,000 repetitions of 4 kips per sq. ft. 
producing a density of 120 p.c.f. and a CBR of 47. Repetition of the test with 
a 2 percent sodium sulfate added to the soil produced a density of 141 and 
increased the CBR to 132. A similar material with a constant moisture con- 
tent of 5.7 percent had a soaked CBR of 90 after 25 blows of the 10-lb. ham- 
mer; after 55 blows the CBR was reduced to 75. However, with 55 blows plus 
7000 repetitions of static load, the CBR was increased to 145. 

Thus, the strength may increase or decrease with increased density de- 
pending upon the structure produced by the method of compaction. 


2. Highway Engr., Physical Research Branch Bureau of Public Roads, U.S. 
Dept. of Commerce, Washington, D.C. 
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George R. Halton, 1 M. ASCE.—The author’s theme is well supported by 
reliable data of the particular type of soil tested. The precautions taken to 
eliminate contributory variables, as exemplified by the use of new material 
for each specimen, are commendable. In many respects the author’s findings 
parallel those of an investigation made in 1943 by the California Division of 
Highways at the Oxnard Flight Strip in Ventura County, California. A few 
facts and observations of the latter investigation, which was conducted under 
the writer’s direction, are given below. 

It will be noted that the Oxnard base material was much coarser than the 
soil from which the author’s data were derived and that the laboratory CBR 
values at densities equivalent to field conditions at Oxnard did not corroborate 
the field performance of the relatively coarse material which failed in spots. 
Otherwise the author’s premise of an inverse strength-density relation above 
certain moisture contents is substantiated by the California study. Both the 
Vicksburg and the Oxnard data show perched saturation and freedom from the 
influence of a high water table. 

The Oxnard Flight Strip was one of many built during the first years of 
World War II under the general supervision of the Public Roads Administra- 
tion with detailed investigation and inspection by state highway departments. 
The material used in its base and subbase was at that time classified as an 
A-2 friable to A-2 plastic gravelly sand. It approximated an SC-GF in con- 
trast to the ML-CL soil at Vicksburg. 

The native subgrade at the Oxnard site was a silty “adobe” clay estimated 
to require a cover of 2.5 to 3.0 feet of 15 CBR subbase and 30 CBR base ma- 
terial. A buried bank of gravelly sand, discovered two miles from the strip 
by the late Col. R. J. Allan?, provided the needed 280,000 cu. yds. of material 
which complied with the CBR and other requirements of the specifications. 

Relative compaction was specified to be not less than 90%, based on the 
California Method? and averaged 96*%, in 118 control tests. The special tests 
of the investigation averaged about 102% The water content was not specified 
but was fairly well regulated, averaging 6.6% in place at the time of sampling 
and 8.8% at optimum. The moisture at the time of rolling was probably one 
or two per cent higher than when sampled. Only 2 of the 118 control tests 
exceeded their respective optimums. 

The above ratio of 2 in 118 is equivalent to the areal proportion of several 
base failures which developed under the trucks hauling bituminous plant- 
mixed leveling course and surfacing. All but a very few of the 45 spots which 
failed had received bituminous surface treatment when distress appeared. 
The defective areas exhibited typical retrogression from an initial slight 
springiness to a final pronounced alligator cracking and rutting. 

During the special investigation to determine the cause of failure, 20 sam- 
ples from failure areas and 6 samples from adjacent areas which showed no 
distress were subjected to tests of field condition and laboratory properties. 
Most of the samples were taken within a few days after damage occurred; a 
few were taken within a few minutes after distress was deliberately induced 
by truck routing. 

Unfortunately neither the equipment nor the precedent for Hveem stabilom- 
eter tests of untreated aggregate or for field CBR tests were locally available 
at that time. The majority of the samples taken in the special investigation 


1. Consulting Engineer, Newark, N.J. 

2. Assoc. M. ASCE; died Feb. 1946; formerly District Materials Engineer, 
Calif. Div. of Highways. 

3. ASTM Procedures for Testing Soils, July 1950, p. 209. 
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TABLE 2 
CCOMP4kISON Ur TESTS AND 
FAILUK™S NON-PAILUKES 
Sample; Sample 
Property Range Avis. | Avy. (ang, Kemarks 


Pass. 3/4" Sieve 5 95 
» Pass.No.59 Sieve 4h 
4 Pass.No.200 Siev 
> Piner than 5 microns 5-9 > 


"> Piner than 1 micron 1-3 


Liquid Limit 14-20 
Plasticity Index 1-6 
PA Class. 

Pileld Conditions 


Density, p.c.f 129-135 131 127 122-129 | (1) (2) 
| G=2.65%0.01 


Nater Content 71-9. 8.3 8.3 5.8-7.9 
Saturation 74-100 74 58-76 (2) 


Laboratory Equivalent 
of rield Condition: 
Nensity, pef fe -132 122-129 
dater Content | 6.1-8.2 
dater Content after 
ye-day soak 7-3-3.3 8.3-9.7 
jain or Loss(=)} 
| #9.5 to 
CHR Before Soakiny | 46-154 > é > 21-196 
CBR After Soaking 89-110 27-92 
Expension % 0,0-9.1 0.1-9.3 
Standard wab. CBK 
Density 122-129 123-126 | (1) (3) 
CBR 65-85 (1) 
| 


In the California Division of Hivhways in 1943, the Casagrande classi- 
fication, the impact compaction of CLK specimens and the adjustment 

of CBR penetration were not standard >rocedure. 

One tailure sanple with apparent density of 140 and saturation of 127> 
was outlawed, 

ASTM Proceaures for Testin,, Soils, July 1950, Pawe 209. 


= 
q a. 
i 
94-96 
3-46 
20-22 
5-7 
2 | 1 2 1-3 
16 14 14 14-16 
Pl. A2 A-2 Fr. | A-2 - (1) 
(1) 
(2) 
(3) 
ay 
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were tested in the laboratory for CBR at the approximate field condition as 
well as at the then standard soaked and unsoaked density achieved under a 
static load of 2000 p.s.i. The results of the special tests are summarized in 
Table 2. 

Infiltration from 0.8 inch of precipitation in the week preceding failure 
was an early cause suspect but was deemed immaterial after field permea- 
bility tests of the bituminous treated upper 3 inches of base material showed 
only 2 x 1074 feet per day. Ground water elevation and capillary action were 
ruled out as causes of failure when the special tests of natural subgrade ma- 
terial showed it to be relatively dry, dense and stable. 

During the repair of the failures the characteristic springiness was ob- 
served in only the upper half of the 2.5 to 3.0- foot combined thickness of base 
and subbase materials. Some of the incipient failures which were not immedi- 
ately repaired subsequently showed no distress when loaded trucks were de - 
liberately routed over them. Elsewhere one or two spot failures were arti- 
ficially induced by deliberate repetition of truck passes. Following the inter- 
ruption due to repairs the paving was completed without the recurrence of old 
failures or the development of new ones. The paved strip was still in good 
condition in 1947, beyond which date the writer has no information. 

The investigation was hampered somewhat by the late Fall paving tempo 
of construction and by the immaturity and divergence of opinions which char- 
acterize problems lacking in precedent. The disparity between field behavior 
and routine or standard CBR values did prompt special tests of laboratory 
CBR at field conditions and in a final field test 3 or 4 weeks after the epidemic 
of spot failures the rear tire of a loaded truck was driven within 8 inches of 
a 6-inch diameter auger hole filled with water in a rubber bag.* The latter 
test produced no observable change in water level. 

The field and laboratory investigations at Oxnard were made by experi- 
enced, expert personnel. Samples from 6-inch diameter auger holes were 
weighed in the field to the nearest 0.01 pound, placed in metal containers with 
sealed covers and checked for weight in the laboratory. Calibrated sand back- 
fill from moisture-proof containers was also weighed to 0.01 pound in the field 
and special depth-density calibrations and techniques assured accurate volume 
determinations at all horizons in the sample holes. At least 50 pounds of ad- 
ditional soil immediately adjacent to the hole was secured from each horizon 
sampled, to provide new material for each laboratory specimen. The entim 
auger portion of each sample was dried at a controlled low temperature and 
the test specimens were moistened for 24 hours before the final addition of 
water required to either match the field condition or straddle optimum, 

It was concluded in 1943 that the Oxnard failures were due to excessive 
saturation in the very dense upper half of base and subbase material and that 
they would have healed or ceased retrogression if left unrepaired for a few 
weeks. The slightly greater plasticity of the samples from failed areas was 
not considered a general cause of failure because 93 of 95 Plasticity Index 
determinations ranged between 0 and 4 because the correlation within that 
narrow spread was inconsistent. The high values of laboratory CBR at equiv- 
alent field condition were attributed to the slow rate of CBR penetration as 
compared to the deformation under a rolling truck; the CBR loading rate al- 
lowed time for the attenuation of excess pore pressures which the truck roll 
compounded into rupture. 

In retrospect these conclusions still seem sound, but the reported differ- 
ences in density and saturation of the failures as against the non-failures 


4. ASTM Procedures for Testing Soils, July 1950. 
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suggest that some unknown factor exaggerated the true differences. The ab- 
normal saturation, the unusual rebound, the difficult laboratory duplication 
of the reported field densities and the loss of water during the soaking period 
all suggest that the auger holes in failed areas were subject to bleeding and 
constriction during the sampling. The inconsistency of computed field den- 
sities on an e-log p chart verifies the suspicions suggested above. 

The qualitative relation of density to pressure is shown in Fig. 6, wherein 
the unit on the log p scale is a theoretical energy ratio rather than a numeri- 
cal quantity of energy per cubic foot (= ft-lbs./ft.3 = unit pressure). The 
virgin slopes were determined from densities achieved in the 5-layer and 
10-layer variations of the California Method. The difference in field densities 
is inconsistent with the uniform, lapped pattern of many roller passes and non- 
uniform, spotty pattern of actual failures. 

The spotty pattern of failed areas suggests that individual loads or groups 
of loads from drier or wetter portions of the pit may have been an antecedent 
cause. Delayed absorption in coarse particles wetted shortly before or dur- 
ing rolling probably caused pore pressures, which were “built in” by the 
roller, to gradually diminish. If the delayed portion of absorption was 0.4% 
by weight = 1.0% by volume, and rolling left 3.0%by volume of air in the voids 
between particles, the total air voids would be 4.0%. If the pressure trapped 
in the 3.0% was 4 tons per square foot, delayed absorption would reduce this 
toward 3 tons per square foot. Compacted soil which experienced such an 
internal relief of pore pressure could then withstand an increment of pressure 
corresponding to 1% of compression, 

Soil without delayed particle absorption would lose trapped pressure, but 
more slowly, via regional diffusion and in proportion with permeabilities and 
lengths of escape paths. Loads applied before such external relief became 
effective could withstand only slight increments of pressure without rupture. 

From the above considerations in conjunction with the author’s data and 
deductions it seems logical that the critical volume of air voids, although of 
the order of 3% on the average, is probably a function of load and loading rate 
as well as a function of the dimensions, type and condition of soil. With be- 
havior subject to so many variables anomalies such as the correlation of field 
CBR with field strength in one soil and the lack of a similar correlation in 
another soil are to be expected. 
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